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Abstract 

General guidelines for constructing a quantum theory of charged-particle 
beam optics starting ab initio from the basic equations of quantum mechanics, 
appropriate to the situation under study. In the context of spin-^ particles, 
these guidelines are used starting with the Dirac equation. The spinor the- 
ory just constructed is used to obtain the transfer maps for normal and skew 
magnetic quadrupoles respectively. As expected the traditional transfer maps 
get modified by the quantum contributions. The classical limit of the quan- 
tum formalism presented, reproduces the well-known Lie algebraic formalism 
of charged-particle beam optics. 
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1 Introduction and Formalism 



Charged-particle beam optics, or the theory of transport of charged-particle beams 
through electromagnetic systems, is traditionally dealt with using classical mechanics. 
This is the case in ion optics, electron microscopy, accelerator physics etc The 
classical treatment of charged-particle beam optics has been extremely successful, in 
the designing and working of numerous optical devices from electron microscopes to 
very large particle accelerators. It is natural to look for a prescription based on the 
quantum theory, since any physical system is quantum at the fundamental level! Such 
a prescription is sure to explain the grand success of the classical theories and may 
also help towards a deeper understanding and designing of certain charged-particle 
beam devices. To date the curiosity to justify the success of the classical theories as 
a limit of a quantum theory has been the main motivation to look for a quantum 
prescription. But, with ever increasing demand for higher luminosities and the need 
for polarized beam accelerators in basic physics, we strongly believe that the quantum 
theories, which up till now were an isolated academic curiosity will have a significant 
role to play in designing and working of such devices. 

It is historically very curious that the, quantum approaches to the charged-particle 
beam optics have been very modest and have a very brief history as pointed out in the 
third volume of the three- volume encyclopaedic text book of Hawkes and Kasper ||. 
In the context of accelerator physics the grand success of the classical theories origi- 
nates from the fact that the de Broglie wavelength of the (high energy) beam particle 
is very small compared to the typical apertures of the cavities in accelerators. This 
and related details have been pointed out in the recent article of Chen ||. 

A beginning of a quantum formalism starting ab initio with the Dirac equation 
was made only recently J7J-0. The formalism of Jagannathanet al was the first one 
to use the Dirac equation to derive the focusing theory of electron lenses, in particular 
for magnetic and electrostatic axially symmetric and quadrupole lenses respectively. 
The formalism of Jagannathan et al further outlined the recipe to obtain a quantum 
theory of aberrations. Details of these and some of the related developments in the 
quantum theory of charged-particle beam optics can be found in the references [|7[-p!4|. 
I shall briefly state the central theme of the quantum formalism with some of the 
more recent results relevant to accelerator optics. 

The starting point to obtain a quantum prescription is to build a theory based on 
the basic equations of quantum mechanics appropriate to the situation under study. 
For situations when either there is no spin or spinor effects are believed to be small 
and ignorable we start with the scalar Klein-Gordon and Schrodinger equations for 
relativistic and nonrelativistic cases respectively. For electrons, protons and other 
spin- 1 particles it is natural to start with the Dirac equation, the equation for spin-| 
particles. In practice we do not have to care about the other (higher spin) equations. 

In many situations the electromagnetic fields are static or can reasonably assumed 
to be static. In many such devices one can further ignore the times of flights which 
are negligible or of not direct interest as the emphasis is more on the profiles of 
the trajectories. The idea is to analyze the evolution of the beam parameters of 
the various individual charged-particle beam optical elements (quadrupoles, bending 
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magnets, • • •) along the optic axis of the system. This in the language of the quantum 
formalism would require to know the evolution of the wavefunction of the beam 
particles as a function of V, the coordinate along the optic axis. Irrespective of the 
starting basic time-dependent equation (Schrodinger, Klein-Gordon, Dirac, • ■ •) the 
first step is to obtain an equation of the form 

d 

ih— ip (x ,y; s) = H (x ,y; s) ip (x ,y; s) , (1) 

where (x, y; s) constitute a curvilinear coordinate system, adapted to the geometry 
system. For systems with straight optic axis, as it is customary we shall choose the 
optic axis to lie along the Z-axis and consequently we have s = z and (x, y\ z) con- 
stitutes a rectilinear coordinate system. Eq. ([[]) is the basic equation in the quantum 
formalism and we call it as the beam-optical equation; Ti and ip as the beam-optical 
Hamiltonian and the beam wavefunction respectively. The second step requires to 
obtain a relationship for any relevant observable {(O) (s)} at the transverse-plane at 
s to the observable {(O) (s in )} at the transverse plane at s in , where s in is some input 
reference point. This is achieved by the integration of the beam-optical equation 
in© 

4>(x,y;s) = U (s,s in )ip(x,y;s in ) , (2) 

which gives the required transfer maps 

(O) (s in ) — (O) (s) = (i(>(x,y,s)\0\1>(x,y;s)) , 

= (^(x,y;s in ) \U ] OU\ tp (x,y; s in )) . (3) 

The two-step algorithm stated above may give an over-simplified picture of the 
quantum formalism than, it actually is. There are several crucial points to be noted. 
The first-step in the algorithm of obtaining the beam-optical equation is not to be 
treated as a mere transformation which eliminates t in preference to a variable s along 
the optic axis. A clever set of transforms are required which not only eliminate the 
variable t in preference to s but also gives us the s-dependent equation which has 
a close physical and mathematical analogy with the original t-dependent equation 
of standard time-dependent quantum mechanics. The imposition of this stringent 
requirement on the construction of the beam-optical equation ensures the execution 
of the second-step of the algorithm. The beam-optical equation is such, that all the 
required rich machinery of quantum mechanics becomes applicable to compute the 
transfer maps characterizing the optical system. This describes the essential scheme 
of obtaining the quantum formalism. Rest is mostly a mathematical detail which is 
built in the powerful algebraic machinery of the algorithm, accompanied with some 
reasonable assumptions and approximations dictated by the physical considerations. 
For instance, a straight optic axis is a reasonable assumption and paraxial approxi- 
mation constitute a justifiable approximation to describe the ideal behaviour. 

Before explicitly looking at the execution of the algorithm leading to the quan- 
tum formalism in the spinor case, we further make note of certain other features. 
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Step-one of the algorithm is achieved by a set of clever transformations and an exact 
expression for the beam-optical Hamiltonian is obtained in the case of Schrodinger, 
Klein-Gordon and Dirac equations respectively, without resorting to any approxima- 
tions! We expect this to be true even in the case of higher-spin equations. The 
approximations are made only at step-two of the algorithm, while integrating the 
beam-optical equation and computing the transfer maps for averages of the beam 
parameters. Existence of approximations in the description of nonlinear behaviour is 
not uncommon and should come as no surprise, afterall the beam optics constitutes a 
nonlinear system. The nature of these approximations can be best summarized in the 
optical terminology as; a systematic procedure of expanding the beam optical Hamil- 
tonian in a power series of |-7Tj_ /j> 1 where po is the design (or average) momentum of 
beam particles moving predominantly along the direction of the optic axis and 7t± 
is the small transverse kinetic momentum. The leading order approximation along 
with |7rj_/po| ^ 1 constitutes the paraxial or ideal behaviour and higher order terms 
in the expansion give rise to the nonlinear or aberrating behaviour. It is seen that 
the paraxial and aberrating behaviour get modified by the quantum contributions 
which are in powers of the de Broglie wavelength (Ao = 2irh/po). Lastly, and impor- 
tantly the question of the classical limit of the quantum formalism; it reproduces the 



well known Lie algebraic formalism |20| of charged-particle beam optics pioneered by 
Dragt et al. 

Let us start with the Dirac equation for a beam made up of particle of charge q, 
mass m and anomalous magnetic moment fi a in a static electromagnetic field with 
potentials (0(r),A(r)) 

R D \$ D ) = E , (4) 

where \ipD) is the time-independent 4-component Dirac spinor, E is the energy of the 
beam particle and the Hamiltonian Ho, including the Pauli term is 

Hd = (3m c 2 + col ■ (-ihV - qA) - p a /5S ■ B , 






= 

*. = (J o ) ' °y = (J o ) ' a * = [ o -i ) ■ (5) 

We multiply H^ (on the left) by a z /c and rearrange the terms to get 

T<d = -poPx®z - qAJ + a z a ± ■ irj_ + (n a /c)f3a z ll ■ B , 

x = { o -r^J' i = ^E~^- (6) 

Eq. (H) is not in the standard form. We need to get Hp into a beam-optical form 
in close analogy with the standard Dirac equation, as required by step-one of the 
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algorithm. Define, 



M = -j={I + X a z ) 



M' 1 = - X a,) 



(7) 



and 
Then 



\^d)-^W) = m\4, d ) . 



d 



ih— = H' , W = MHdM" 1 = - Po f3 + s + o. 



(8) 
(9) 



with the matrix elements of 8 and O given by 



En 

£12 
£22 



- q A z i - (ji a /2c) {(z + r 1 ) <r±-B ± + (z- r 1 ) a,B z } 



So, = 



- q A z t - (ji a /2c) {(z + r 1 ) o- _l • b ± - (e - r 1 ) ^} , (10) 



and 



On 

012 



0. 



22 



21 



• ± ■ *± - (//./2c) {i (f - r 1 ) - B^) 

- (e + r 1 ) s,n}" , 
r 1 [«r ± • tt ± + (/i a /2 C ) {i (e - r 1 ) (5,^ - s,^) 

+ (e+r 1 )sji}; ■ 



(ii) 



These transformations gives us the beam-optical Hamiltonian as required by step- 
one of the algorithm. In the Dirac theory, the lower spinor components are much 
smaller than the upper spinor components. This is also true for the beam-optical 
equation derived above. The analogy between the standard Dirac equation and the 
beam-optical equation, thus constructed can be summarized in the following table 



Standard Dirac Equation 

m c 2 p + £ D + 6 D 
m c 2 

Positive Energy 
Nonrelativistic, \tt\ -C itiqc 
Non relativistic Motion 
+ Relativistic Corrections 



Beam Optical Form 

- Po f3 + £ + 6 
-Po = 

Forward Propagation 
Paraxial Beam, |7Tj_| <C hk 
Paraxial Behavior 

+ Aberration Corrections 



From the above table it is clear that there is a very close analogy between the de- 
rived beam-optical form and the standard Dirac equation. Having established this 
analogy we can now execute step-two of the algorithm. As is well known the Foldy- 
Wouthuysen machinery |2T] enables one to take the proper nonrelativistic limit of 
the Dirac equation and provides a systematic procedure for obtaining the relativistic 
corrections to a desired degree of accuracy in powers of 1/moc 2 . With the analogy in 
the table above, we are now in position to adopt the Foldy-Wouthuysen machinery to 
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the beam-optical situation. The procedure of the Foldy-Wouthuysen-like machinery 
gives us the paraxial behaviour and a systematic procedure to compute aberrations 
to all orders in powers of 1/po- To leading order we get 



|V/)=exp(— - (3d) i> 



(i) 



(12) 



Then 



ih 



d_ 

dz 



exp 



(i) 



i 



V 2p 
ifrexp 



(3d ft'exp — (3d 



1 



£W=£ -^(3Q 2 + 
2p 



1 



2p 

(3d) ^- I exp f — /30 



dz 



0(D 



2 P{ 



0,8 



2 Po 







(13) 



The lower pair of components of are almost vanishing compared to the upper 

pair and the odd part of liS x > is negligible compared to its even part we can effectively 
introduce a Pauli-like two-component spinor formalism based on the representation. 

Calling as H 



m— 

dz 



n 



n 



1> 



qA z + — Trj 
2p 



{(q + e)B z S z + 7 eB ± -S ± } 



With 7T i 



^ 2 _L ^ 2 



e = 2m fi a /h, 'j = E/m c 2 , S = -her . (14) 



Up to now, all the observables, the field components, time etc., have been defined 
in the laboratory frame. The covariant description of the spin of the Dirac particle has 
simple operator representation in the rest frame of the particle. In the analysis of the 



spin-dynamics in accelerator physics it is customary to use such a representation [18 



So we define the following transform which takes us from the beam-optical form to 
the accelerator optical form 



V>) = exp j^- (Tf x a y - fr y o- x yf • 



(15) 



Details of the accelerator optical transform are found elsewhere [p~3j , [L6|, [17]. Up to 
the paraxial approximation the accelerator-optical Hamiltonian is 



ih- 

dz 



■0 



(A) 



H(A) 



(A) 
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V 2 Po J Po 

with O s = — tqB + efBu + jB ± )\ . (16) 

7mo 1 v n 

Only at this stage we need to know the specific geometry of the magnetic lens, i.e., the 
form of the fields characterizing the lens, for computing the transfer maps. We shall 
briefly consider the normal and skew quadrupoles respectively in the next sections 
and see how the traditional transfer maps get modified by the quantum contributions. 



2 Example-I: Normal Magnetic Quadrupole 

An ideal normal magnetic quadrupole of length I is characterized by the field 

B = {-Gy,-Gx,0), (17) 
corresponding to the vector potential 



A=(o,0,iG(x 2 -y 2 )) 



(18) 



situated in the transverse-planes z = z- ia and z = z ont + i, with G constant inside the 
lens and zero outside. The accelerator-optical Hamiltonian is 



H 



H(z) 



f = -Po + ^P± , for z< z in and z > z ont , 
H L (z) = -po + ^pl - \qG (x 2 - y 2 ) + 3f (ya x + xa y ) , 

2 



(19) 



for z- m < z < z, 



out ) 



with rj = (q + je)G£h/2pl . 



The subscripts F and L indicate, the field-free and the lens region respectively. 

Best way to compute U is via the interaction picture, used in the Lie algebraic 
formulation [21J of classical beam optics. Using the transfer operator thus derived (de- 
tails found elsewhere fL3| , |17| ) we get the transfer maps for averages with the subscripts 
in and out standing for (z in ) and (z out ) respectively 



( ( x ) ) 




/ 


( (x) ) 


( 


(Px)/P0 






(Px}/P0 














(y) 






(y) 




\ (Pv)/Po ) 


out 


V 


V (Pv)/Po ) 


\ 



\ 



^sinh (VJCi 



'K I 
cos (VkA-1 
Kl 
sin (VKi 



fa) 



n 



f cosh( v / ^Q ^sinh(V^) 

~Ksmh(^fKl) cosh{y/K£) 



V 



in 



cos(v / ^^) ^ sin( v / ^^) 
-VKsm(y/K£) cos(\fK£) J 
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M< 

> 



/ 1 Az< \ 

1 > 

1 Az< 

V o o o i ) 



T 



M>M q M< 



(20) 



For spin, the transfer map reads 



(Sy) 



yl out 



out 



('S'x)in 

Ann 

(Sy) in 

Attt] 

Wz)in 



sin(V^) 
iff 



'cos(VFl) - I s 



/cosh( v / ^^) - l\ n , /cos( v / ^Q - l\ n , I 

+ I I (p x S x ) in + ( (pA)in • (21) 



Kipo 



K£p 



Thus, we get the a fully quantum mechanical derivation for the traditional transfer 
maps (transfer matrices) ||]. In addition we also get the spinor contributions, the 
Stern-Gerlach kicks as they are called. In recent years there has been a campaign to 
make a spin-splitter ||2!| device to produce polarized beams using the Stern-Gerlach 
kicks. The spinor theory of charged-particle beam optics, in principle supports the 
spin-splitter devices. 



3 Example-II: Skew Magnetic Quadrupole 

For a skew-magnetic-quadrupole lens the field is given by 

B = (-G 3 y,G s x,0), (22) 
corresponding to the vector potential 

A = (0, 0, —G a xy) . (23) 
The accelerator-optical Hamiltonian is 

{H F = ~Po + 2^P± , for z < z in and z > z out , 
H L (z) = -p + ±pl + \qG s xy - (xa y - ya x ) , (24) 
for z- m < z < z out , with n s = (q + >ye)G s £h/2pl . 

and the corresponding transfer operator is 

U i t L^Z on t y 2qn) 
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I 1 Vs 



S + C~ A 5- (C+ - 2) \ 
-p x + -^r~ n Vx + -7=- n PoV H — — Py cr y 



s- 



nP0X + 



-Px 



C± = cos(JK s £) ± cosh( JKJ), S ± = sm(JK s £) ± wahUK, 



(25) 



Then we get the transfer maps 

/ (x) \ (( (x) \ 



(Px)/Po 

(y) 

(Pv)/Po 



\ 



sq 



' out 



M. 



sq 



s+ 



(Px)/Po 

(y) 

(Py)/Po 

c~ 



+ V_s 

2 



/ -C-(a y ) + (C+-2)(a x ) \\ 
K„£ 

-S+{a y )+S-{a x } 
VK~ s t 

-(C+-2)K)+C-(^) 
KJ 



J 



V 



-5-(a i/ )+5+(a 3: ) 



JK S S- C H 



fK~ s S+ C- 



//in 

T sq = M>M sq M< . (26) 



V -v^5 + a- -vk s s~ C + J 

Again we get a quantum derivation of the traditional transfer maps (3| with some 
spinor contributions. Due to the algebraic machinery of the quantum formalism it is 
straightforward to see that the transfer maps of the skew quadrupole can be obtained 
by rotation of the transfer maps for the normal quadrupole by 7r/4 along the optic 
axis. For spin, the transfer map reads 



{Sy)out ~ (Sy)m + 



(4), 



(Sm)i 



Ao 



An 



2ttt] s 



5^ 



S~ 



s- 



5 M 



+ 



s- 



+ 



s- 



fK. 



5 H 



j (yS g ) in + 
) (a^in + 

;J (^)i„ + 

(xS x )i n I 

^ (2/5r)in - ( 
(ySy)in + I 



c- 



K s £p Q J 
( C + -2 
\ K s £p 
' C+-2 N 

( C~ > 

K s £p 0/ 

(—) 

' C+-2 

c- 



K. 



(PxS z )m 
J (Py&}taj , 
m 

\Pa; 'S'a; ) in 

(Py*S , x)in|' 
J (Px^)in 
(pA>inH- (27) 



9 



4 Concluding Remarks and Directions for Future 
Research 



To summarize, an algorithm is presented for constructing a quantum theory of charged- 
particle beam optics, starting ah initio from the basic equations of quantum mechanics 
appropriate to the situation. This algorithm is used to construct a spinor theory of 
accelerator optics starting ab initio from the Dirac equation, taking into account the 
anomalous magnetic moment. The formalism is demonstrated by working out the 
examples of the normal and skew magnetic quadrupoles respectively and as expected 
there are small quantum contributions, which are proportional to the powers of the 
de Broglie wavelength. It is further shown that the quantum formalism presented, 
in the classical limit reproduces the Lie algebraic formalism [50] of charged-particle 
beam optics. 

The present algorithm is suited for constructing a quantum theory of charged- 
particle beam optics at the single-particle level. The next logical step would be 
to extend it to accommodate the realistic multiparticle dynamics. We feel that 
clues to such a formalism can come from the experience gained from the so-called 



quantum- like models p3| . These phenomenological models have been extensively de- 
veloped in recent years to explain the classical collective behavior of a charged-particle 
beam, by associating with the classical beam a quantum- like wavefunction obeying a 
Schrodinger-^&e equation with the role of h being played by the beam emittance e. 

One practical application of the quantum formalism would be to get a deeper 
understanding of the polarized beams. A proposal to produce polarized beams using 
the proposed spin-splitter devices based on the classical Stern-Gerlach kicks has been 
presented recently pZf . 

Lastly it is speculated that the quantum theory of charged-particle beam optics 
will be able to resolve the choice of the position operator in the Dirac theory and the 
related question of the form of the force experienced by a charged-particle in external 



electromagnetic fields p4j| . This will be possible provided one can do an extremely 
high precision experiment to detect the small differences arising in the transfer maps 
from the different choices of the position operators. These differences shall be very 
small, i.e., proportional to powers of the de Broglie wavelength. It is the extremely 
small magnitude of these minute differences which makes the exercise so challenging 
and speculative! 
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